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ABSTRACT 


New  temperature-compensated  high  piezoelectric  coupling  crystals  are 
needed  for  acoustic  surface-wave  matched  filters,  encoders,  decoders,  and 
related  signal  processing  devices.  Elastic  coefficients  ,  together  with  their 
pressure  and  temperature  derivatives,  are  also  important  in  materials  science 
and  in  geophysics.  Oxides  show  several  correlations  between  structure  type 
and  elastic  behavior.  Stiffness  coefficients  of  silicate  minerals  are  generally 
larger  in  the  direction  of  Si-0  bonding.  A  simple  mechanical  analogy  in  which 
atomic  bonds  are  simulated  by  springs  connected  in  series  and  parallel  is  used 
to  estimate  the  size  and  anisotropy  of  the  elastic  moduli.  Magnitudes  of  the 
temperature  and  pressure  derivatives  of  the  stiffness  coefficients  are  also 
predicted  by  the  model.  Crystallochemical  factors  conducive  to  anomalous 
elastic  behavior  are  examined,  with  special  attention  to  materials  which  beco&e 
Stiff er  with  increasing  temperature.  Internal  rotational  motions  involving 
bending  of  the  atomic  bonds  appear  particularly  important  in  this  class  of 
material*. 
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INTRODUCTION 

Recent  development  of  acousto-optic  and  surface-wave  devices  has  re¬ 
newed  interest  in  the  elastic  properties  of  oxides.  New  temperature-compensated 
materials  with  large  piezoelectric  moduli  are  especially  needed.  In  a  delay- 
line  medium,  an  electric  signal  is  converted  to  an  acoustic  wave  which  travels 
at  a  slower  velocity  and  is  subsequently  re-converted  to  an  electric  signal. 

Elastic  coefficients  govern  the  velocity  of  the  acoustic  wave,  and  hence  the 
transit  time. 

Elastic  properties  of  oxides  are  also  important  in  the  geosciences, 

•Inca  most  minerals  in  the  earth's  crust  are  oxides.  Stiffness  moduli,  bear 
special  relevance  to  our  understanding  of  rock  mechanics  and  seismic  wave  velocities 

Prom  previous  work  it  appears  that  packing  density  is  the  primary  var¬ 
iable  affecting  the  elastic  moduli  of  oxide  compounds.  Birch  (1961a,  1961b)  shoved 
that  most  common  minerals  have  about  the  sarnie  mean  atomic  weight  (molecular 
weight  divided  by  the  number  of  atoms  in  the  chemical  formula),  and  that 
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longitudinal  sound  velocity  is  roughly  proportional  to  density.  Shear 
velocities  also  increase  with  rock  density  (Simmons,  1964) .  Anderson  and 
Nafe  (1965)  plotted  hulk  modulus  (reciprocal  volume  compressibility) 
against  volume  per  ion  pair,  demonstrating  that  most  oxide  data  follow  the 
same  relation.  Bulk  modulus  is  inversely  propprtional  to  volume,  regardless 
of  whether  the  volume  change  is  caused  by  juicrsure,  compositional  variation, 
temperature,  porosity,  or  phase  changes.  As  might  be  expected,  bulk  modulus 
increases  with  density  because  short-range  repulsive  forces  make  it  increasingly 
difficult  to  compress  the  solid  as  the  atoms  move  closer  together. 

STRUCTURE-STIFFNESS  CORRELATION 


All  solids  change  shape  under  forces.  Under  small  stress,  the  strain  e 
is  related  to  the 'stress  0  by  Hooke’s  law,  (a)  »  (c)(e).  The  elastic  stiffness 
coefficients  (c)  constitute  a  fourth  rank  tensor  in  which  the  number  of 
independent  coefficients  depends  on  symmetry  (Nye,  1957).  In  contracted  matrix 
notation  0^  *  cijej  w^,ere  i»  j  ■  1»2,3  refer  to  longitudinal  stresses  and 
strains  along  axes  X^,  Xj  and  respectively.  For  shearing  notions  about 

each  of  the  axes,  i,  j  »  4,5,6.  The  discussion  which  follows  will  be  concerned 
with  the  relative  values  of  c^„  Cj2  and  the  coefficients  relating  change 
in  length  along  a  principal  direction  to  a  parallel  component  of  Stress. 

The  stiffness  coefficients  for  silicate  minerals  presented  in  Table  1 
show  a  correlation  between  elastic  anisotropy  and  structure  type.  All 
silicates  contain  (SiO^)  tetrahedra  and,  depending  on  how  the  tetrahedra  are 
linked  together,  con  he  classified  as  framework,  layer,  chain  or  ring  silicates. 

In  framework  silicates  such  as  quarts  and  silica  glass,  the  tetrahedra 
form  three-dimensional  networks.  Since  the  bonding  is  nearly  isotropic,  there 
is  no  cleavage  and  littld  anisotropy  in  hardness  or  elasticity.  Compare  the 
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Table  1 

v 

'  Longitudinal  Stiffness  Constants  Expressed  in  Megabars  (I1EARM0N,  1966) 


Framework  Silicates 


Cl-quartz 

Si02 

C11,B  c22  °  0,9 '  c33  “  2,2 

silica  glass 

Si02 

C11  ’  c22  -  c33  -  °'8 

Single  Chain  Silicates  (pyroxenes) 

aegerine 

NaFeSi206 

C11  ’  2'S-  c22  '  1'8-  ‘  2'3 

augite 

(CaMgFe)Si03 

C11  *  2,8>  c22  “  1'5'  c33  "  2,2 

diops ide 

CaMgSigOg 

ca  -  2.0,  c22  -  l.S,  c33  -  2.4 

Double  Chain  Silicates  (amphiboies) 

hornblende  (Ca  .tta.K)  2.3  WtJc>A1)  5(sl  ,41)  8<>22  WU)  2 

eu  .'l.2;  Cj2  i  1.8,  eJ3V  M 

Ring  Silicates 

beryl  BejAljSigC^g  C11  "  c22  "  3.1  c33  -  2.8 

tourmaline  (Na,Ca) <Li,Mg, Ai) 3{Al,FofMn) 6(0U> /t  (803)  3S*6°18 

C11  ’  '»  "  2-2'  c33  ‘  1-6 

taygr  Silicates 


biotite 

KCMg.Fe)  3(A1S  iplQ)  (Oil)  % 

C11  “ 

c22  “  1,$» 

c33 

muscovite 

K  Al2(A\Si3Oio)<OH)2 

cll  * 

c22  *  l*8* 

C33 

phlogopitc 

K  Hg3(AlSi3010K0U)2 

Cll* 

c22  * 

C33 

-  0.5 


-  0.6 


-  0. 
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longitudinal  elastic  moduli  given  in  Table  1.  When  corrected  for  density 
(p,  «  2.65  gms/c.c.  for  quartz  and  2.2  gms/c.c.  for  silica  glass)  the 
stiffness  constants  are  nearly  identical  for  the  two  forms  of  SiC^,  sub¬ 
stantiating  the  relation  between  bulk  modulus  and  volume  (Anderson  and 
Nafe,  1965). 

The  influence  of  crystal  structure  becomes  more  obvious  in  the  chain 
silicates.  Pyroxenes  contain  SiO^  single  chains,  and  amphiboles  Si^O^  double 
chains  ns  shown  in  Fig.  1.  Elastic  coefficients  in  Table  1  are  referred  to 
the  measurement  directions  denoted  by  arrows  in  Fig.  1.  Bonding  is  stronger 
along  the  chain  direction  giving  rise  to  pronounced  cleavage.  We  also  expect 
the  crystal  to  be  stiffer  in  this  direction,  resulting  in  larger  moduli. 
Experiment  confirms  this  suggestion;  the  stiffness  parallel  to  the  chain  (c^j) 
is  the  largest  in  pyroxenes  and  amphiboles.  We  also  not©  that  has  increased 
considerably  in  Hornblende,  possibly  because  of  the  increase  of  chain  width 
in  this  direction. 


Beryl  and  tourmaline,  two  ring  silicates,  show  a  similar  correlation 
between  stiffness  and  structure.  Both  contain  rings  as  illustrated 

schematically  in  Fig.  1.  We  expect  strong  bonding  and  greater  stiffness  in 


the  plane  of  the  ring,  hence  should  be  smaller  than  and  as  observed. 
Beryl  is  not  very  anisotropic  because  of  the  strong  Bc-0  and  Al-0  bonds 


connecting  the  rings. 


When  Si60lg 


rings  adjoin  one  another,  the  tetrahedral  layer  found  in 


micas  Is  formed.  The  cleavage  and  stiffness  anisotropy  become  very  obvious  in 


lays"  silicates,  where  c^  and  are  three  times  larger  than  c^j.  This  is 


the  maximum  elastic  anisotropy  observed,  for  reasons  that  are  explained  later. 


Fig.  1.  Arrangements  of  SiCV  tetrahedra  in  silicates. 

(a)  Single  chain  silicates,  (b)  Double  chain 
silicates,  (c)  King  silicates,  (d)  Layer 
silicates. 
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MECHANICAL  ANALOG 

-  la  the  lattice  tiu-ory  of  elastic  coefficients,  stiffness  coefficients 
are  related  to  atomic  force  constants  hy  determining  the  energy  associated 
with  various  strain  components.  The  calculation  is  cumbersome  for  a  mon¬ 
atomic  simple  cable  lattice  Chattel,  1953),  and  would  be  overwhelming  for 
most  mineral  structures.  To  avoid  mathematical  complexity  and  gain  further 
physical  insight  regarding  the  causes  of  elastic  anisotropy,  we  make  use  of  a 
simple  mechanical  system. 

The  analogy  used  to  describe  elastic  anisotropy  is  one  in  which  two 
mechanical  springs  represent  atomic  bonds  with  force  constants  k  and  K.  Tables 
of  force  constants  (Wilson,  Decius  and  Cross,  1955)  derived  from  infrared 

4 

vibrational  spectra  show  that  typical  values  are  within  an  order  of  magnitude 

of  10  dynfes/cm  (1  md/A).  Stretching  foree  constants  are  plotr* !  as  a 

function  of  interatomic  ui stance  in  Fig.  2.  As  might  be  <«pccted,  the  force 

constant  decreases  with  increasing  bond  length,  though  it  is  rather  surprising 

that  so  many  different  typos  fall  along  the  same  lino.  In  general  the  strongest  bonds 

have  the  largest  force  constants  because  deep  potential  wens  have  larger  second 
derivatives  when  well-shapes  arc  similar.  Tims,  for  example,  it  is  found  that  the 
stretching  force  constants  for  C-C,  OC,  and  OC  are  about  5,  10  and  16  ed/A, 

respectively,  For  bending  ooticuc,  the  force  csr.s tarts  li**’  the  range 
0.1  -  0.6  md/A.  Bending  force  constants  are  such  smaller  than  stretching  force 
constants  because  the  repulsive  overlap  potential  dominates  when  a  toes  awe 
toward  ono  another,  increasing  the  force  constant  for  stretching  action. 

To  explain  the  clastic  properties  of  solids  containing  both  strong  and 
weak  bonds,  consider  the  spring  system  illustrated  in  Fig.  3.  When  strong 
and  weak  springs  ave  connected  in  series,  most  of  the  elastic  energy  i® 
stored  in  the  weak  springs,  while  in  the  parallel  connection  the  strong  spring 


STRETCHING  FORCE  CONSTANT  (md/R) 
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Fig,  2.  Relationship  between  force  constant  and  interatomic  distance.  Stretching 
force  constants  are  from  Wilson  Decius  and  Cross  (1955),  and  interatomic 
distances  from  Sutton  (1958) , 


1.0  1.2  8.4  1.6  1.8  2.0  2.2  24 

BOND  LENGTH  (X) 


.  fig.  3.  Series  and  parallel  connections  of  springs  used 
to  represent  atomic  bonds.  Strong  bond3  have 
large  force  constants  K,  while  weak  bonds  are 
easily  stretched  and  have  smaller  constants  k. 


contains  most  of  the  energy.  Let  K  and  k  be  the  force  constants  of  two  bonds 
arranged  in  series  and  parallel  positions,  as  shown.  This  is  a  schematic 
representation  cf  the  bonding  in  mica.  In  muscovite  Si-0  and  K-0  bonds  are  in 
series  for  tensile  stresses  applied  perpendicular  to  the  sheet  and  in  parallel 
whan  the  applied  forces  lie  in  the  plane.  In  pyroxenes,  the  parallel  connection 

applies  to  measurements  along  the  silicate  chains,  and  series  connections  to  the 
two  perpendicular  direction. 

Analyzing  the  series  arrangement  for  an  applied  tensile  force  F  gives' 

F  »  aA  *  n  Km.  +  n  ku,  «*  t/'u  where  0  is  the  stress  acting  on  a  surface 

q  S  R  S  K  5  8  S  o 

of  cross-sectional  area  A  containing  n  chains,  u  /l  is  the  resulting  strain 
in  the  series  (s)  connection  whose  overall  stiffness  is  cg.  u^  and  u^  are  the 
displacements  of  the  springs  with  force  constants  K  and  k.  An  identical  force 
applied  to  the  parallel  arrangement  gives  an  analogous  expression.  F  ■  crAp  » 
n  Ku  +  n  ku.  *  c  A  u  ll  .  For  the  series  connection  both  springs  experience 

pK  pkpppp 

the  same  force  so  that  their  restoring  forces  are  equal,  uRK  «  u^k.  The  total 
displacement  ug  -  ufe  +  uR,  giving  u^  *  ug/(l  +  k/K)  and  u^  -  ug/(l  +  K/k),  and 


Vs 

“V 


v  k+K  ' 


(1) 


In  solving  the  parallel  chain  system,  it  is  obvious  that  the  displacements  of 
the  different  springs  are  equal  and  that  up  «  2u^  *  2^.  Substitution  in  the 
force  equation  gives  the  elastic  constant 


n  t 
Q  ■  ■■  J-JP. 
P 


s 


(2) 


Note  that  cg  and  c^  are  unequal,  even  when  all  the  springs  ere  identical  (k  •  K). 
The  elastic  coefficients  depend  on  bond  lengths  through  l  and  on  the  number  of 


chains  per  unit  area  in  different  directions 
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To  determine  the  effects  of  strong  and  weak  bonding  on  the  elastic  con¬ 
stants,  assume  that  the  geometric  factors  are  about  equal  so  that  n  l  /A 
n8£8/\'  giving  the  ratio  .Cg/c^  »  4kX/(k+K)  .  In  Pig.  4,  the  quantity  cg/Cp  is 
plotted  as  a  function  of  K/k  to  illustrate  the  effects  of  mixed  bonding.  When 
K/k  «  1,  all  bonds  have  the  same  force  constant,  and  the  elastic  constants  are 

of  course  the  same  for  series  and  parallel  connection,  so  that  c  /c  ■1. 

8  P 

At  the  other  extreme  cg/c^  G  as  K/k  00  ,  but  the  approach  to  zero  is  very 
slow.  For  K/k  **  2,c  «  0.9c  and  even  for  K/k  «  10,  c  «*  0.3c  .  Since  force 
constants  for  various  chemical  bonds  are  all  within  an  order  of  magnitude  of 
one  another,  the  expected  maximum  elastic  anisotropy  is  about  3:1,  as  observed 
in  muoco/ito  (Table  1).  In  all  cases  c  >  c  as  observed  experimentally. 

r  ° 

It  would  be  useful  to  predict  the  magnitude  of  the  observed  stiffness 
coefficients  (Table  1)  as  well  as  their  anisotropy.  For  chains  identical  in 
size,  number  and  force  constants,  equations  (1)  and  (2)  reduce  to 


r  nf, 

Q-  m  c  *  k  » 

p  a  A 


n,  lt  and  A  ran  be  evaluated  from  the  crystal  structure,  but  k  cannot.  Typical 
values  for  k are  0.1  to  1.0  mJllidynes/A  for  bending  force  constants  and  1-10 
cd/X  for  stretching  constants.  Both  type°  of  deformation  come  into  play 
in  minerala.  When  an  SIO^  tetrahedron  is  stressed,  fot  instance  both  stretching 
and  bending  will  take  place.  Si-0  force  constants  (Matossi,  1949)  arc  fairly 

typical  with  stretching  constants  4-5  tnd/A  and  bending  constnntsO.6  to  0.9ad/X, 
wh^lo  those  involving  A1  are  10-20%  smaller  (Hidalgo  and  Sorratosa,  1956)*  A 
‘careful  analysis  would  be  required  to  dot  ermine  the  correct  force  constant  to 


Stiffness  anisotropy  for  series  and  parallel 
connections  plotted  as  a  function  of  spring 
constant  ratio  K/k.  The  series  connection  ii 
far  more  pliant  when  K>>k. 


use  in  (3) >  but  in  any  case  all  the  atomic  force  constants  have  not  been 
determined  from  spectroscopic  data.  When  an  average  value  k  VL  md/A  is  sub¬ 
stituted  in  (3) ,  elastic  stiffness  cc  jfficients  of  the  right  magnitude  are 

obtained.  Taking  1  'v  3 A,  —  ^  1/2,^  ^QlA  ^  gives  c  ^3 /x.  10^*  -(!Xn-?s.  »  3  mega- 

cm 

bars,  comparable  to  the  experimental  values  (Table  1) .  * 


PRESSURE  DEPENDENCE  OF  THE  ELASTIC  STIFFNESS 
*  The  pressure  derivatives  of  the  elastic  coefficients  of  minerals  determine 
changes  in  seismic  wave  velocities  deep  within  the  earth,  and  are  strong 
indicators  of  the  onset  of  phase  transformation.  Elastic  stiffness  coefficients 
ahd  their  initial  pressure  derivatives  for  fo  ir  oxides  are  listed  in  Table  2» 
Using  the  spring  model  just  described,  three  observations  are  to  be 
rationalized: 

(1)  The  pressure  derivatives  aro  all  about  one  to  ten  megabar/mogabor. 
(dimensionless) . 


(2)  Largo  stiffnesses  usually  show  greater  pressure  derivatives  than 


small  ones:  if  >  .then  *gp^“  >  gp — ♦ 

(3)  Pressure  derivatives  of  the  stiffnesses  arc  positive  in  dense-packed 


structures  but  in  open  structures  are  occasionally  negative,  iiuarti  ona 
beryl  each  have  one  negative  derivative  but  the  close-packed  corundum  and 
forsterite  structures  show  none. 


To  estimate  the  pressure  dependence  of  the  elastic  stiffness  we  again 

2  ‘ 
make  the  approximation  that  n/A  •v  1/f.  ,  then  c  ■»  k/1,  and 


Ac/A3P  -  (1  }t)  (Ak /£&)  -  ( kAt2)  (At/AP) ,  Assuming  an  isotropic  solid,  AP  v.  (At)  (c ID 


'able  2 
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Comparison  of  plastic,  stiffness  and  their  initial  pressure  derivatives 
for  four  oxides*  Adiabatic  stiffness  c^j  are  expressed  in  megabars, 
and  the  pressure  derivatives  3c^j/3P  are  dimension] ess. 


Bex\vl  Quartz  Corundum  Forsterlte 


ij 

cu 

cu 

aV3p 

cu  . 

Scjj/SP 

CU 

aCij/3F 

11 

3.09 

4.5 

0.87 

3.3 

4.98 

6.2 

3,29 

8.3 

22 

It 

t< 

II 

It 

•1 

it 

2.01 

5.9 

33 

2.83 

3.4 

1.06 

10.8 

5.02 

5.0 

2.36 

6.2 

44 

0.66 

-0.2 

0.58 

2.7 

1.47 

2.2  • 

0  .67 

2.1 

55 

II 

u 

II 

»i 

«i 

ii 

0.81 

1.7 

66 

0.90 

0.3 

0.40 

-2.7 

1.68 

1.5 

0  .81 

2.3 

*  12 

1.29 

3.9 

0.07 

8.7 

1.63 

3.3 

0  ,07 

4.3 

13 

1.19 

3.3 

0.12 

6.0 

1.17 

3.7 

0  *07 

4.2 

23 

M 

it 

« 

M 

•  * 

it 

« 

0  .07 

3.5 

14 

■m:  • 

-0.18 

1.9 

-0.23 

0.1 

' 

.«»• 

Reference  Yoon,  1971  HcSkirain,  Gieske  and  Grahrna  and 

Andrea tch6  Barsch,  Barsch,  1969 

Thurston,  ■  1968 

1965 
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A  rough  value  for  Ak/A&  can  be  obtained  by  examining  how  k  varies  with  ^ 

interatomic  distance.  Short  strong  bonds  have  larger  stiffnesses  than  long 

bonds.  The  bond  stiffness  for  Si-0  is  about  10%  larger  than  that  for  Al-0, 

and  the  bond  length  is  about  10%  shorter.  Therefore  Ak/Ai  is  roughly-2rad/A2. 

Stretching  force  constants  are  plotted  as  a  function  of  interatomic  distance 

in  Pig.  2.  For  distances  greater  than  1.5A,  the  relation  between  k  and  <£  is 

°2 

nearly  linear  with  a  slope  of  about  -2  md/A  . 

Substituting  this  value  in  (4)  along  with  c  **  3  x  10^2  dyne/cm2,  k  *  1  md/A 
and  Z  a  3&  gives  Ac /AP  *v  10 (dimensionless) ,  the  right  order  of  magnitude. 

To  explain  the  second  observation  consider  the  anisotropic  structure  in 

l 

Pig.  5a,  The  structure  contains  tightly  bonded  atoms  in  the  X^  direction  and 
very  loose  bonding  along  X2.  From  arguments  previously  presented 
Now  consider  their  pressure  derivatives.  From  eq.  (3)  the  change  in  stiffness 
with  pressure  is  related  to  the  change  in  the  number  of  chains  per  unit  area 
n/A,  their  repeat  distance  and  bond  stiffness  k: 

Ac  „  A(n/A)  .  nk  Aft  .  nft  Ak  fSV 

Under  pressure  the  structure  will  compress  mainly  along  because  of  the  weak 
bonding  in  that  direction,  giving  the  exaggerated  deformation  in  Fig.  5b, 

For  the  direction  there  will  be  little  change  in  and  so  that 
<ACjj,/AP)<vfcJki  A(n/A)/AP.  The  number  of  chains  per  unit  area  increases  rapidly 
With  pressure  because  of  the  big  reduction  in  fcj*  decreasing  A  and  increasing 
n/A.  Therefore  increases  rapidly  with  pressure. 
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For  direction  X2,  there  is  little  change  in  n/A  with  pressure  because 

hardly  changes.  Therefore  Ac  /AP  z  (n/A)  [k,  (A£  /AP)  +  t0  Ak  /AP) .  The 
1  22  1  2  2 

length  and  spring  constant  are  inversely  related  to  one  another  so  that  the 

increase  in  spring  sonstant  is  partially  offset  by  change  in  length.  Hence  C22 

will  not  increase  rapidly  with  pressure. 

The  third  observation  regarding  the  pressure  dependence  is  the  occurrence 
of  negative  derivatives  in  open  structures  like  beryl  ami  .quarts.  Whan  a 
dose-packed  structure  is  compressed,  the  atoms  move  clou«*r  together  but  this 
need  not  be  true  in  an  open  structure  where  rotations  can  take  place.  To 
determine  the  effect  on  the  elastic  constants,  consider  equation  .(5)  describing 
the  pressure  dependence  of  the  stiffness.  The  new  feature  here  is  the  pressure 
dependence  of  the  stiffness  k.  If  we  are  considering  the  stiffness  along  X^, 
for  example,  k^  may  decrease  with  P  because  at  high  pressure  a  stress  along 
Xj  produces  a  bonding  rather  than  a  stretching  motion,  the  stiffness  co¬ 
efficients  for  bending  are  considerably  smaller  than  for  stretching.  Thus 
rotation  can  lead  to  negative  pressure  dependence  of  shearing  stiffness  co- 
■  efficients. 

Temperature  dependence  of  the  elastic  stiffness 

Host  materials  soften as  they  become  warmer*  so  (hot  the  elastic  stiffness 
decreases  as  temperature  Increases.  Some  representative  values  of  dc/d T for 
oxides  are  given  in  Table  3.  tae  fractional  change  in  stiffness  (1/c)  (dc/dT) 
is  of  the  order  -2  x  10  per  °C.  The  variation  of  the  elastic  coefficients 
with  temperature  is  of  some  geophysical  interest  because  of  the  increase  in 
temperature  with  depth*  Stiffness  increases  with  pressure  so  that  acoustic  waves 
increase  speed  with  depth*  but  the  increase  is  partly  offset  by  the  temperature 
effect,  the  temperature  dependence  of  the  elastic  constants  is  also  important 


Table  3.  Temperature  coefficients  of  elastic  stiffness  for  several  oxides 

-4 

T  ■  (l/c)(dc/dT)  and  is  expressed  in  units  of  10  per  °C. 
c 


MgO 

Tc^  *  -2.3 

Xc44  -  -1.0 

smo3 

Tcu  *  -2-6 

Tc. .  »  -1.1 
44 

osio2 

Tcu  -  -0.5 

Tc33  “ 

Tc44  '  -!•« 

Tc66  "  +1*6 
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in  acoustic  delay  lines  and  in  piezoelectric  oscillators  and  filters.  Variations 

of  the  delay  time  or  resonant  frequencies  are  undesireable  in  these  devices. 

The  mechanical  spring  model  just  discussed  also  provides  an  explanation  of 

the  temperature  dependence  of  the  elastic  stiffness.  To  estimate  the  size  of 

dc/dT,  we  assume  an  isotropic  model  of  identical  springs,  for  which  c  ■  n  Z  k/A  ~ 

k ft,  as  shown  previously,  k  is  the  force  constant  of  the  spring  and  Z  the  length. 

The  temperature  derivative  is 

dc  d_ 
dT  13  dT 


<!> 


l 

I 


dk 

dT 


k,  dl 
'rif. 


The  second  term  is  directly  proportional  to  the  linear  thermal  expansion  coefficient 
o  «  (1 ft)  (d£/dT) .  The  first  term  depends  on  the  change  in  force  cons  taut  with 
temperature.  Chemical  bonds  grow  weaker  with  increasing  length  and  the  force  con¬ 
stant  becomes  smaller.  If  we  assume  that  the  change  in  k  depends  primarily  on  lt 
then  (dk/dT)  «  (dk/dv) (d£/dT),  and  the  variation  stiffness  with  temperature  becomes s. 


dk  k 

d l  “Z 


a 


6}  «  9  t  <» 

Aft  shown  in  the  previous  section,  dk/df  -  -2  od/A  ,  k, - 1  ed/A,  f v  3A,  and  0.3x10 
dynes/eo  .  Thermal  expansion  coefficients  for  oxides  are  about  10  per  aC, 


Substituting  these  values  gives  a  fractional  change  (l/c)(dc/dT)  «►  10  /°C, 
which  is  the  right  order  of  magnitude.  This  value  is  small,  but  is  several 
times  larger  than  the  thermal  expansion  coefficient  a.  This  is  important 
because  the  thermal  expansion  coefficients  arc  crucial  in  delay  line  and 
frequency-standard  devices,  thus  in  searching  for  zero-temperature  coefficients, 
materials  with  unusual  elastic  properties  are  needed.  In  particular,  it  is 
important  to  find  materials  with  positive  temperature  coefficients  U/c) (dc/dT). 
In  the  next  section  ve  consider  the  circumstances  under  which  this  sight  occur . 
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Temperature  -  Compensated  Materials 

The  transit  time  t  in  a  delay-line  medium  is  given  by  the  transit  distance 
&  divided  by  the  acoustic  velocity  /c/p  .  c  is  the  appropriate  stiffness  co¬ 
efficient  and  p  the  density.  Temperature  variations  alter  the  transit  time, 
causing  problems.  The  temperature  derivative  (l/t)(dt/dT)  *  (l/£)(d£/dT)  - 
(1/2  c)(dc/dT)  +  (1/2  p)(dp/dT).  The  first  term  on  the  right  is  the  linear 
thermal  expansion  coefficient  a,  measured  parallel  to  the  transit  direction. 

The  third  term  ia  also  related  to  a.  For  an  isotropic  material,  (l/p)(dp/dT)  «  -3a 
Thus  the  temperature  drift  of  the  delay  time  is  governed  by  the  thermal  expansion 
coefficients  and  by  the  change  in  Stiffness  with  temperatures 
Tt  -  (1/tKdt/dT)  £  -l/2(a  +  Tc). 

T  ie  (l/c)(dc/dT) ,  the  fractional  change  in  stiffness  coefficient  with  temperature 

v 

and  is  approximately  ~2xXQ  ^/*C,  as  discussed  in  the  previous  section.  Thermal 

expansion  coefficients  arc  usually  positive  and  somewhat  snailer,  about  +10  /°C. 

The  temperature  variation  of  the  transit  time  is  therofore  detemiaud  by  X  , 

c 

and  is  about  10  Observed  vatuea  of  Tt  reported  for  surface  wave  devices 

are  slightly  smaller,  about  30-90  ppm/°C  (Carr,  1972). 

Temperature-compensation  requires  that  Tfc,  the  temperature  coefficient  of 
the  transit  time,  be  tero.  For  most  materials,  Tfc  is  positive  because  elastic 
stiff ness  decreases  with  Increasing  temperature,  slowing  down  tbs  acoustic  waves 
and  lengthening  the  transit  time.  Density  also  decrease  with  temperature,  speeding 
up  the  waves  and  shortening  the  transit  time,  but  this  offset  is  generally  smaller 
then  that  due  to  the  elastic  constants* 

Materials  with  unusual  values  of  T&  are  therefore  of  special  interest  for 
temperature-compensated  devices.  Quarts  is  such  material  since  T  for  sheer  about 

V  .■ 

the  t-exls  is  positive  (Table  3).  The  structures  of  a  and  0-  quarts  (Fig,  6.) 

* 

provide  an  understanding  of  this  behavior.  The  a -quart*  structure  ia  a  partially- 
collapsed  derivative  of  0-quart* .  At  higher  temperatures  the  SiO^  tetrehedra 


2C 


Fig.  6.  Arrangement  of  silicon  atoms  in  quartz.  Oxygens  are  located  approximately 
midway  between  neighboring  silicon  atoms. 


Of  —QUARTZ 


-“QUARTZ 
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rotate  to.  an  open,  fully-expanded  structure,  undergoing  a  diaplacive  phase 
transition  at  3?3°G  to  tue  grquartz  structure.  Quarts  shows  a  rather  high  rate  of 
thermal  expansion  at  room  temperature  due  to  the  rotation  of  the  fcetrahedra. 

Thermal  expansion  coefficients  become  even  larger  just  below,  the  trans- 
formaxion  temperature  and  then  decrease  abruptly  to  near  zero  whan  the 'fully-  ... 
expanded  structure  is  achieved.  Quantitative  calculations  (Taylor,  1972) , 
show  that  the  thermal  expansion  of  quartz  is  due  chiefly  to  tetrhhedral  rotations, 
with  only  a  minor  contribution  from  thermal  motion.  The  anomalous  value' of  Tctg' 
is'  related  to  the  rotation  effect,  e^g. relates  a  shearing  stress  about  the 
Z'-axis  to  a  shearing  strain  about  the  ..same  axis  -  the  same  type  of  action  involved 
in  the  rotation  of  tetrahedra  in  transforming  from  a  to  $  quartz.  At  room 
temperature  the  Si  -  0  -  Si  bonds  are  quite  bent,  (Fig.  7  a),  but  they  straighL’eri' 
with  increasing  temper ture  (Fig.  7 b)  as  the  structure  grows  closer  to  (Quartz. 
Consider  the  effect  of  mechanical  stress  on  such  arrangements,  When  the,  bent 
structure,  is  stretched,  the  fcetrahedra  rotate  as  well  as  deform  (Fig.  7c) ,  resulting 
in  a  sizeable  strain.  Thus  the  stiffness  Is  rather  small  at  room  temperature. 

When  the  high-temperature  structure  is  stressed  (Fig.  7d) ,  only  deformation  occurs 
since  the  structure  is  too  straight  to  permit  rotation.  As  a  result,  there  is 
less  strain  and  the  stiffness  increases  with  temperature.  A  positive  value  of 
Tr  can  be  relat-  i  to  rotational  effects  in  this  way. 

In  an  earlier  section,  the  temperature  coefficient  of  the  elastic  stiffness  c 
was  described  in  terms  of  a  mechanical  spring  model.  There  it  was  shown  that 
Tc“  (l/c)('’c/dT)  -  (l/£c)  (dk/dT)  -  (k/^2c)  (d£/dT) . 

I-  this  equation,  l  is  the  length  of  the  spring  and  k  is  the  spring  constant. 

For  most  mater?  ,1s,  T  is  determined  by  the  first  tarm  on  the  right,  resulting 

c 

from  the  change  in  sprin^  constant  with  temperature.  The  second  term  arises  from 

thermal  expansion,  and  is  generally  smaller.  The  spring  constant  usually  decreases 

with  temperature  because  of  thermal  vibrations,  so  that  dk/dT  is  usually  negative 

and  so  is  Tc  . 

\ 


Fig*  ?,  i^dels  indicating  fcfe' origin  ®f- positive  t^  coeffieiencs  in  open  "crystal 
£ true  Lures,  (a>  Shows  three  atoms  connected  by  bent  bonds.  At  higher  temperatures 
(k).r;£})'S''  «uni  stvMghr.ens Avhcm  stressed,  the  bent '  bonds  both  lengthen  and  rotate, 
L&&.  shown  in  (t)  ,  '.WheiiFtbe  straight  .bonds  are  stressed,  only  lengthening  takes 
$lace  (<?)  >-  (c)-shews  a  larger- dveirall  change  in  length  than  (d) ,  hence  the  lov- 
teracratura  ottuatn-e  is  more  compliant.  The  bent  bonds  in  a-quartsr  are  shown 
in  {a>  ,  *  A;  s  h^hjL'tiig  stress  about  X^  pr  duces  tension  along  the  bond  direction, 
resulting ‘in. band  rotation ; as  veil  as  extension.  For  the  high  temperature 
^qusrta'strdOiurQ  (f),'!np  bending  occurs  ih  this  prelection. 


HIGH  T 


O  silicon 
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This  need  not  be  so  when  rotations  occur  because  rotational  spring  constants 
are  smaller  then  stretching  opring  constants.  For  Si-0  bonds,  the  stretching 

o  » 

constant  is  about  5  md/A  while  bending  spring  constants  are  less  than  1  md/A. 

In  the  kinked  bonds  (Fig.  6a)  both  stretching  and  bending  occurs  under  stress 
(Fig.  7c).  The  effective  force  constant  is  smaller -than  in  the  high- temperature 
case  (Figs.  7b  and  7d)  where  only  stretching  takes  place.  In  this  situation 
dk/dT  is  positive,  as  is  T  ,.  because  the  force  constant  changes  character  with 
increasing  temperature.  ' 

Paratellurite,  one  of  the  polymorphs  of  tellurium  dioxide,  exhibits  unusual 
elastic  behavior  similar  to  quartz.  Temperature-compensated  cuts  of  para- 
tellurite  have  been  used  in  surface  wave  devices,  although  the  piezoelectric 
coupling  parameters  are  small  (Carr,  1972).  In  TeC^  the  shearing  stiffness 
cg  **  1/2  (c^  -  c12)  has  a  positive  temperature  coefficient.  The  crystal 
structures  of  the  paratellurite  and  rutile  polymorphs  of  TCO^  are  illustrated 
in  Fig.  8.  Paratellurite  is  a  partially-collapsed  derivative  of  the  rutile 
form.  When  paratellurite  is  heated,  the  bent  configuration  straightens, 
similar  to  the  changes  observed  in  quartz.  Thus  the  shearing  spring  constant 
increases  with  temperature  making  dc/dT  positive. 

In  searching  for  new  materials  with  positive  temperature  coefficients,  we 
look  to  open  structures  in  which  sizeable  rotations  can  occur.  Coordination 
numbers  are  usually  small  in  open  structures,  so  that  oxygens  are  bonded  to 
only  1,  2,  or  3  neighbors.  Close-packed  structures  are  generally  not  conducive 
to  rotational  movements,  although  there  «ie  exceptions  to  this.  SrTiO^  shows 
a  beautiful  rotational  transition  near  110°K  which  is  accompanied  by  anomalous 
elastic  behavior.  The  Ti-O-Ti  bonds  are  kinked  in  the  low- temperature  structure. 
Leucite,  eristobalite,  anorthite,  and  other  silicate  minerals  show  ’puckering*  • 
transitions  similar  to  quartz  but  unfortunately  these  materials  are  not  piezoelectric. 
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Phase  transitions  are  a  second  clue  to  the  occurrence  of  positive  T 

V 

values,  since  rotations  sometimes  lead  to  a  change  in  structure.  ’Poled *  ferro- 
electrics  such  as  LilaO^  or  ferroelastics  are  therefore  of  Interest.  To  pre¬ 
vent  domain  formation,  the  transition  should  not  be  near  room  temperature.  A 
material  such  as  Bi^WO^.  (T,  ■  950°C)  might  be  of  interest.  The  elastic  constants 
of  bismuth  tungstats  have  not  been  determined.  BigWOg  is  polar  at  room  temp¬ 
erature  but  not  ferroelectric.  Like  lithium  niobate,  it  is  a  'frozen*  ferro¬ 
electric. 

Anomalous  elastic  properties  are  also  observed  in  BajNaNb^O^^,  known 
colloquially  as  bananas.  The  compliance  coefficients  and  s55  decrease 
with  increasing  temperature  (Yamada,  Iwasaki  and  Niiseki,  1970).  As  with  SIQ2 
and  TeO-,  the  unusual  behavior  is  associated  with  a  phase  transition  to  higher 
symmetry.  Barium  sodium  niobate  undergoes  a  displaciye  phase  transition  at 
300* C  which  is  primarily  elastic  in  character.  The  x  and  y  exes  of  the  room- 
temperature  orthorhombic  structure  are  rotated  by  45*  from  those  of  the  higb- 
temperature  tetragonal  phase,  Ba2Nahb^0^g  exhibits  large  piezoelectric 
coupling  constants  in  certain  directions. 
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